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1. intrqducao 


A necessidade de se preservar a vida util dos reservato 
rios de agua, tern levado o honiein a procurar ineios niais capazes de cons 
tatar e control ar os fatores que possani anieagar a longevidade desses 
corpos d'agua. Dentre estes fatores, uni dos niais importante sao as par 
tvculas solidas em suspensao, transportadas por interniedio de sous tri^ 
butarios e acuinuladas no corpo principal da represa. 

Conforiiie Ritchie et al (1976), as tocnicas utilizadas 
normalniente para a medida dos sedinientos eni suspensao na agua, envoi vcm 
a medida da sua turbidez por nieios oticos ou gravimetricos, que deter 
minam sua concentragao de sedimentos. Essas iiiedidas de acordo com es_ 
tes mesmos autores, sao de utilidade limitada, desde que el as represen 
tarn pontos de amostragem e desta forma nao perniitem uma visao global 
de toda a area da represa. 

Di versos estudos tern demons trado haver uma alta corre 
lagao entre a quantidade de sedimentos em suspensao e os nTveis de re 
flectancia da radiagao solar, iiiedidos por interniedio de sistemas de sen 
sorianiento remote e, que estes nTveis permitiriam uma estimativa quan 
titativa da concentragao destes sedinientos. 

A vantagem de se utilizer as niedidas estimadas por inter 
medio desses sensores a bordo de satelites, como os da serie LANDSAT 
ou aeronaves, seria a visao sinotica fornecida pelas iniagens obtidas 
por esses sistemas. 0 ideal, na realidade, seria a conjugagao dos dois 
sistemas de medidas, o convencional e o obtido por sonsoriaiiiento remo 
to, de modo a medir precisamente a quantidade de sedimentos e a sua 
distribuigao no reservatorio de agua. 
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Este trabalho tern conio objetivo mostrar a utilidade da 
analise automatica de imagens MSS do satelite LANDSAT, para caracte 
rlzar a distribuigao das partTculas solidas era suspensao na superfTcie 
da represa de Tres Marias e estimar a sua concentragao na superfTcie 
da agua. 


Os resuUados aqui apresentados foram obtidos numestudo 
que esta sendo desenvolvido pelo Instituto de Pesquisas Espaciais jun 
tamente com a Companhia de Desenvolviinento do Vale do Sao Francisco, na 
area da represa de Tres Marias e seu sistema tributario, e representani 
apenas parte do mesmo, 

2. MATERIAL E METODOS 


A represa de Tres Marias esta local izada na parte cen 

tral do Estado de Minas Gerais (Figura 1), apresentando uma area inun 
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dada de 1102 km . Ela esta situada numa regiao onde predominam latO£ 
solos vernielho escuro e verraelho amarelo distroficos (EPAMI6, 1978), 
formados sobre rochas de Pleistoce No - Holoceno e Pre-Cambriano (CPRM, 
1977). 

As fitas magneticas, conipatTveis com coniputador - CCT 
(Computer Compatible Tape), correspondentes a orbita 164, ponto 25, da 
passagem de 5 de abril de 1978, foram utilizadas para compartimentar 
a represa em diferentes areas, de acordo com as suas respostas espe£ 
trais, apresentadas em termos de nTveis de cinza, atraves de urn proces 
so automatico supervisionado e interativo, no sistema IMAGE-100(I. 100). 

Os dados de reflectancia gravados nas fitas CCT, foram 
armazenadas na memoria do sistema I-lOO, com nTvel de resolugao igual 
a 256, por intermedio do programa Entrada e Saida do VTdeo ("VTdeo 
I/O"). Simultaneamente estes valores sao apresentados pelo sistema na 
tela do Video colorido, formando uma imagem compos ta. 
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A area desta imagem composta, correspondente a represa 
de Tres Marias, situada entre as coordenadas da imagem des 

x-j * 1 390 J 3s"i 

Xg 3186 yg e 1422 

foi ampliada por intermedio do Programa Carregador de CCT para as eSM 
las: 


horizontal * 1:356,907 
vertical ** 1:401.380 

e gravada no Disco 6, no arquivo denominado REPRES, 478. 

Em seguida. fez-se a delimitaclo da area referente, ex 
clusivamente a lamina d 'igua da represa. Para isto foratn selecionadas 
amostras de treinamento, por intermedio do cursor, sobre a area da re 
presa definida claramente na cena do canal MSS 7. Com escolha dessas 
areas e a aplicagao do Programa de Celula Onica ("Single Cell"), forani 
definidos os limites nnninio e maximo de variagao da rcsposta espectral 
da agua neste canal, os quais sao apresentados a seguir: 

Limite inferior = 1 
Limite superior = 28 

Feita a delimitacao, gravou-se esta area no arquivo de 
nominado REPRES. LIM, do disco 0 (zero) 

De modo a permitir a visual izagao do espago espectral 
da assinatura, a ser obtido a partir da aplicagao do programa de Maxi^ 
ma Verossimilhanga (MAXVER) (Figura 2 e 3), foi necessario comprimir 
os valores de nivel de cinza de 256 para 64, Assim, utilizou-se o pro 
grama "COMPRES", tomando-se a area da represa definida no arquivo 
REPRES. LIM, como area de interesse (Figura 4), 
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pressao a ser utilizado no arquivo repres, 478 (Di£ 
CO 6), usando comp area de treinamento o tema Repres, 
Lim (Disco 0 e 6). 



Etn seguida usou-se o progratra "CONTRAST STRECH , nos ca 
nals MSS 4 e 5. de modo a se obter uma melhor visual izagao da comparti 
mentagao da represa, para a escolha das areas de treinamento. A sele 
gao das olto areas de treinamento, fol realizada atraves do resultado 
da apl Icagao deste programa no canal MSS 5, pelo fato do mesmo posslbj^ 
lltar uma melhor visual Izagao da compartimentagao da represa, em dife 
rentes nTveis de cinza (Figura 5). 


Areas de treinamento selecionadas a partir 
dos diferentes nfvels de cinza existentes 
na lamina d'agua. 


^ 
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Para cada area de treinamento foi realizada uma classi 
ficagao do tipo Celula Om'ca ("SINGLE-CEU”), considerando-se apenas 
os canals MSS 4 e 5. Obtendo-se desta forma a distribui^ao de frequen 
da dos nTveis de cinza de cada area de treinamento. A apresentagao des 
ses histogranias de frequencia permite verificar a forma da distribui 
gao dos dados (Figura 6, 7, 8, 9, 10, 11, IP, 13). No caso desta nao 
ser Gaussiana, caracterTstica basica para a aplicagao do ProgramadeMa 
xima Verosimilhanga, e necessario que se proceda a urn novo treinamento 
e classi ficagao ou que se faga urn processo de refinamento de assinatui 
ra, ate que se satisfaga esta condigao. 

Apos este procedimento foi aplicado o programa MAXVER, 
usando, como dados basicos, as 8 areas do treinamento previamente esco 
Ihidas, 


0 programa MAXVER, destina-se a classificagao ponto a 
ponto de imageni multi ‘•espectrais* obtidas em geral de satelites. 0 me 
todo de se utilizer este programa pode ser dividido em duas fases dis 
tintas: a) de treinamento, na qua! o usuario escolhe e determina os pa 
rametros espectrais; e b) de classificagao propriamente dita. 

Como resultado, pode-se ter estatTsticas para cada cla£ 
se de treinamento, alem da propria imagem classificada ponto a ponto 
(Velasco et al, 1978). 

Os resultados da classificagao feita atraves do si sterna 
MAXVER, foram apresentados na forma de uma imagem tematica colorida, 
na tela do video col ori do, e na forma grafica, por intermedio da impre^ 
sora de linhas. 
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Fig. 7 - Distribuigao de frequencia de nTveis de cinza na classe 2. 
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de frequencia de nTveis de cinza na classe 3 
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Fig. 11 - Distribuigao de frequencia de nTveis de cinza na classe 6 
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Fig. 12 - Distribuigao de frequencta de niveis de cinza na classG 7 
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Ao resultado desta classificaijao, foraw associad'ts os va 
lores da profundidade de Secchi» oxpressa em metros, que foram obtidas 
durante o trabaVho de campo realizado nesta represa, no perTodo de 31 
de margo a 19 de abril de 1978. Estes valores foram medidos na superfT 
cie, em 24 pontos de amostragem, conforme indrcado na Figura 14. 

Estes pontos de amostragem foram definidos a parti r de 
uma classificagao tematica feita em uma fita CCT, da passagem de 06/ 
02/77. Para a escolha dos pontos de amostragem, foi levado em conside 
ragao os diferentes nTveis de cinza, encontrados na superfTcie da re 
presa, e os principais rios abastecedores da mesma. Em cada area da re 
presa que houvesse urn nivel de cinza representative e urn rio abastece 
dor de importancia, foi colocado urn ponto de amostragem, de modo a per 
mitir que se fizesse uma estimativa da concentragao de partTculas soli^ 
das em suspensao, em cada classe mapeada. 

3. RESULTADOS E DISCUSSAO 


A partir da analise automatica das fitas CCT, foi obti 
da uma classificagao tematica da represa, mostrando a compartimentagao 
da mesma em oito classes diferenciadas em tennos de nTveis de cinza. U 
to foi apresentado no video colorido do sistema I-lOO e, por interme 
dio de urn "Slide”, obteve-se urn negative fotografico colorido, cuja co 
pia em papel fotografico e apresentada na Figura 15. 

Pela impressora de linhas obteve-se urn mapa tematico nas 
escalas aproximadas: Horizontal - 4:157.583 VeJ'tical - 1 :156.176, 
que e apresentado em escala reduzida na Figura 16. 

As classes representadas no mapa e na classificagao te 
matica, apresentada no video colorido, foram definidas pelo programa 
MAXVER, a partir dos valores de media e matriz de covariancia, calcuU 
dos com os dados de resposta espectral encontrados nas areas de treina 
mento, cujos os resultados sao apresentados na Tabola 1. 
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As classes foram ordenadas dc acordo com a intensidade 
da resposta espectral, atrlbuindo-se para cteitos visuais, coloracao 
azul para a classe com a resposta espectral mais haixa e vcrmelho para 
a resposta espectral mais alta. As classes intormcdiarias receberam 
uma combina^ao de cores, indo do azul ate o vermelho, passando pclo ver 
de (Figura 17). 0 mapa obedece a esta mesma ordena^ao. 



Classif icagao temati a da reprcsa de Ire. 
Marias, atraves do si sterna MAXVf.R, com a 
gcragao e ordenagao de ( ores do acordo i om 
os nTveis de cuiza. 
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A label a 2 Mostra como foram gerados as cores para cada 
uma das classes. 


TABELA 2 

GERACAO DE cores para as 8 CLASSES DETER MINADAS 
ATRAVES DO PROGRAMA MAXVER 


ntveis 

RGB 

para 

classe 

1 

> 

0 

0 

255 

NIVEIS 

RGB 

para 

classe 

2 

> 

0 

110 

255 

NTVEIS 

RGB 

para 

classe 

3 

> 

0 

255 

255 

NTVEIS 

RGB 

para 

classe 

4 

> 

0 

255 


NTVEIS 

RGB 

para 

classe 

5 

> 

0 

255 

110 

NIVEIS 

R G B 

para 

classe 

6 

> 

llj! 

1 255 1 

NTVEIS 

RGB 

para 

classe 

7 

> 

255 0 

255 

ntveis 

RGB 

para 

classe 

8 

> 

23/3/55 


R - Vermel ho 
6 - verde 
B - azul 

Atraves da analise do mapa, pode-se observar que as re£ 
postas espectrais mais baixas local izarani-se proxiinas a saTda de agua 
e as respostas espectrais mais altas, local izaram-se nas entradas de 
aguas principalmente aquelas alimentadas pelos rios Paraopeba, Sao Fran 
cisco e Borrachudo. 

Sobre este mapa foram local izados os 24 pontos de aiiios_ 
tragem, onde foram realizadas as medidas de profundidade Secchi , permi_ 
tindo-se fazer uma relagao das respostas com os va lores de profund i_ 
dade obtidos, os quais sao apresentados na label a 3. 
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TABULA 3 

RELACAO DOS PONTOS DE AMOSTRAGEM COM AS RESPECT! VAS 


MEDIDAS DE PROFUNOIDADE SECCHI 











































A Tabela 4 apresenta a rela^ao existente entre as respos 
tas espectrais e as medldas de profundidade Secchi . 

TABELA 4 

MEDIAS DAS RESPOSTAS ESPECTRAIS NOS CANAIS MSS 4 e 5 


E MEDIDAS OE PROFUNDIDADE SECCHI EM METROS 


CLASSES 


MEDIAS 

MSS4 

MSS5 

16.81 

4.72 

20.31 

9.24 

22.11 

9.20 

26.46 

16.47 

30.47 

23.69 

1 

33.64 

33.17 

35.51 

34.02 

36.41 

51.66 


PROFUNDIDADE EM M 


4.00 


1.75 a 3.00 


2.25 


1,75 a 2.00 


1.25 a 1.50 


1.00 a 1.25 


1.25 


0.25 a 0.80 


























Como pode ser observado na label a 5, as classes 2,3 e 
4 poderlam ser agrupadas em uma unica classe, o niesnK) acontecendo com 
as classes 6 e 7. 

TABELA 5 

RELACAO das classes, MEDIAS ESPECTRAIS, PROFUNDIDAOES 


SECCHI E PONTOS DE AMOSTRAGEM 


CLASSES 


MSS4 


16.81 


2 20.31 9.24 


22.11 


26.46 16.47 


30.47 23.69 


33.64 33.17 


35.51 34.02 


36.41 51.66 


PROFUND IDADE EM M 

PONTOS DE AMOSTRAGEM 

4.00 

2A - 2C 

1.75 a 3.00 

2B,1B,4B i 

2.25 

3A 

1.75 a 2.00 

3B, 3C, 5A 

1.25 a 1.50 

4A, 4C, 1C, 8A 

1,00 a 1.25 

7B, 8B 

1.25 

7A 

0.25 a 0.80 

6A, 6B, 6C, 7C 
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0 agrupamento das classes 2,3 e 4 apresentaria as seguin 
tes caracterTsticas: 

a) resposta espectral media no canal MSS4, variando de 20.31 a 
26.46; 

b) resposta espectral media no canal MSS5, variando de 9.24 a 
16.47; 

c) profundidade Secchi, variando de 1.75 a 3.00 m. 

Ja no agrupamento das classes 6 e 7, as caracterTsticas 
seriam as seguintes: 

a) respostas espectral media no canal MSS4, variando de 33.64 a 
35.51; 

b) resposta espectral media no canal MSS5, variando de 33.17 a 
34.02; 

c) profundidade Secchi, variam do de 1.00 a 1.25 ni. 

Este fa to pode ser explicado pel a local iza?ao das areas 
de treinamento, selecionadas para representar essas classes, ocorrereni 
em zonas de transigao sendo, portanto, as mesnias pouco diferenciadas 
entre si, deficultando a definigao de seus limites. 

A Figura 18 mostra o resultado das unices das classes, 
feita no si sterna I -100. 
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Deve-se salientar que os pontos de aniostragem forain e£ 
tabelecidos a partir de unia arialise automatica prelimlnar', de uma fita 
CCT, do mes de fevereiro de 1977, e que os mesnios foram mantidos para 
0 perTodo em que foi realizado o trabalho de campo deste estudo. !sto 
se deve ao fato de que a referida fita era a unica idformagao do perTo 
do chuvoso, existente ate a epooa do trabalho de campo, uma vez que a 
imagem utilizada neste estudo, so foi obtida pelo sensor MSS do LANDSAT 
por ocasiao do mencionado trabalho. 

Devido a este fato, os pontos de amostragem podem nao 
estar localizados com precisao, de modo a serem perfeitamente repre^ 
sentativos das classes por eles caracterizados. 

Com foi mencionado anteriormente, deve-se observar que 
este trabalho apresenta apenas os resultados parciais de uma pesquisa 
mais ampla, que esta sendo desenvolvida pelo INPE e CODEVASF, envolvein 
do dois perTodos sazonais em tres anos diferentes. 

0 ponto de amostragem mais problematico foi oponto lA, 
que apresentou resultados anormais. A explicagao para este acontecimen 
to podera ser dada atraves de sua observagao em outros anos, no mesmo 
perTodo sazonal. 

Urn fator causou problemas na interpretagao automatica 
das fitas CCT, foi a cobertura de nuvens sobre alguns bragos da repre 
sa. 


0 resultado da comparagao desses dados mostrou que os 
niveis de cinza, que representam as diferentes reflectancias da agua, 
sao inversamente proporcionais a profundidade Secchi, isto e, quanto 
maior for a reflectancia, menor e a profundidade Secchi, Como a pro^ 
fundidade Secchi (SD) esta correlacionada ao coeficiente de atenuagao 
da agua (a) pel a relagao: 
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a 

e este fa tor, por sua vez, esta correlacionado com a quantidade de pa£ 
tTculas em suspensao na agua, conclul-se que, conhecendo>se a relagao 
matematica entre a profundidade Secchi e a quantidade de particulas so 
lldas em suspensao, poder-se-ia estabelecer valores quantitativos da 
dispersao de sedimentos em suspensao, para as diferentes classes de nT 
vels de cinza. 

Para exemplificar tal afirmagao, usaremos dados da Ta^ 
bela 5 para obtermos o coeficiente de atenuagao das classes 1 e 8. 

1 - Exemplo da Classe 1; 

Reflectancia - MSS4: 16.81 
MSS5; 4.72 

Ponto de amostragem - 2A 

Profundidade Secchi - 4.00 m 

SDriil = 0,57 

a SO 4,0 


2 - Exemplo da Classe 8: 

Reflectancia - MSS4 - 36.41 
MSS5 - 51.66 


Ponto de Amostragem - 6C 
Profundidade Secchi - 0.25 m 



2,3 ^ 

SO 0.25 


a 


a = 


9,2 
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Como pode ser observado no exemplo acima, as reflectan 
cias da classe 1, o coeficiente de absorgao do urn dos pontos de amos 
tragem pertencentes a esta classe e a sua profundidade Secchi sao bem 
menores do que os da classe 8. Portanto, podemos dizor que a quantida 
de de sedimentos em suspensao na classe 1, ponto 2A, e bem menondoque 
na classe 8, ponto 6C. 

Conclui-se, portanto, que, sabendo-se a profundidade de 
Secchi, podemos calcular o coeficiente de atenua^ao da agua. Associan 
do-se este as diferentes classes de nTveis de cinza, pode-se ter o 
comportamento da dispersao de sedimentos em susnensao na superficie 
d'agua, alem de se obter, para o mesmo, valores quantitativos nas difo^ 
rentes classes de nTveis de cinza, 

4, CONCLUSAO 


Este trabalho mostrou que as respostas espoctrais de cor 
pos d'agua, gravadas nas fita CCT, obtidas a pnrtir do dados inedidos 
pelo sensor MSS do satelite LANDSAT, podem ser usadas no atual estagio 
de pesquisa do INPE, para fazer o mapa de contorno da distribuigao das 
partTculas solidas em suspensao, atraves de sua analise automatica no 
si sterna I -100. 


A utilizagao rotineira desses dados, permitiria a cons 
tatagao rapida e eficiente de algum problema em grande escala, que vie£ 
se a acontecer na represa, relacionada principalmente a distribuigao 
das partTculas solidas em suspensao. 

Alem disso, esse mapa poderia ser utilizado para estab£ 
lecer pontos de amostragem mais representatives, capazes de permitirem 
0 controle da qualidade da agua de reservatorios semelhantes a este. 

A conjugagao desses mapas com urn si sterna de amostragem 
eficiente e contTnuo, permitiria obter dados quantitativos e consequeji 
temente urn resultado mais efetivo para este tipo de estudo. 


I 
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